Strongly Correlated Topological Superconductors and Topological Phase Transitions via Green's 
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We propose several topological order parameters expressed in terms of Green's function at zero frequency 
for topological superconductors, which generalizes the previous work for interacting insulators. The coefficient 
in topological field theory is expressed in terms of zero frequency Green's function. We also study topological 
phase transition beyond noninteracting limit in this zero frequency Green's function approach. 
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I. INTRODUCTION 

Recently topological phases have been among the central 
topics in condensed matter physics^^. These phases of mat- 
ter are interesting in that their bulk is gapped but their sur- 
face is robustly gapless. Among these topological phases 
are topological insulators, including quantum Hall insulators^, 
and the recently discovered time reversal invariant topolog- 
ical insulators including the quantum spin Hall insulators'^*^ 
and its three-dimensional(3d) generalization'^. On the other 
hand, many superconductors have gapped fermionic spectrum 
in the bulk, therefore, they can also be topologically classi- 
fied. For free fermion systems, the topological classification 
of insulators and superconductors are well established by pre- 
vious works^i^. The central tools in these topological clas- 
sifications are various band topological invariants ^-^'^2^^, 
starting from the classical Thouless-Kohmoto- Nightingale- 
den Nijs(TKNN) invariant"'. In the presence of interactions, 
topological states of matter can be generally defined in terms 
of the topological response functions within topological field 
theory''. Recently, there are great interests in the field of in- 
teracting topological insulators Among the most urgent 
problems in the field of interacting topological insulators is to 
formulate topological invariants/topological order parameters 
to distinguish different topological phases. To this end, topo- 
logical order parameters expressed in terms of Green's func- 
tion for interacting topological insulators were proposed 
There are much recent interest in the Green's function ap- 
proach to interacting topological insulators^""'''*. More re- 
cently, it was proposed that topological order parameters can 
be expressed in terms of Green's functions at zero frequency 
without losing any information^^-^*, which greatly simplifies 
numerical calculations^**-^^'^^** and analytical calculation. 

The theory of interacting topological superconductors are 
still less investigated compared to topological insulators. In 
j^gf 39,40 jj. jj^g been proposed that gravitational topological 
field theory can be formulated to describe interacting 3d time 
reversal symmetric topological superconductors how- 
ever, there is no explicit formula for the 9 angle appearing 
in the 6 term (even for noninteracting systems). Therefore, 
it is generally a difficult task to determine the topological 
class of a given interacting superconductor. Recently, there 
are several studies on interacting superconductors in 2d^^ 
in which some judiciously chosen symmetries are imposed 



instead of the more conventional symmetries such as time re- 
versal symmetry. In the present paper, we will focus only on 
the more conventional cases, namely superconductors with- 
out any additional symmetry and superconductors with time 
reversal symmetry, since these cases are closer to real mate- 
rial problems in experiments. By extending the previously 
mentioned Green's function approach to superconductors, we 
will obtain several simple yet general topological order pa- 
rameters for interacting superconductors. These topological 
order parameters are expressed in terms of zero frequency 
Green's function, which are easy to compute by various nu- 
merical methods and analytical calculations. We also present 
simple examples in which interaction effect induce nontriv- 
ial topological phase transitions, which are nevertheless well 
described within our approach. 

The rest of this paper is organized as follows. In SecJIIjwe 
present topological order parameter for 2d chiral supercon- 
ductors expressed in terms of zero frequency Green's func- 
tion, and then in Sec|III] we study a simple but nontrivial 
ansatz to illustrate the application of our new formula in 
strongly correlated systems. In SecJIVI we present topologi- 
cal order parameter for 3d time reversal invariant topological 
superconductors. Then in Sec|V]we obtain the gravitational 9 
angle in terms of Green's function, and show that it is exactly 
given by the winding number of zero frequency Green's func- 
tion proposed in SeclIVI In Sec lVIl we discuss generalization 
of zero frequency Green's function approach to gapless sys- 
tems. We then make several concluding remarks in Sec JVIIl 
including finite size calculations and the general recipe of the 
zero frequency Green's function approach to topological in- 
sulators/supercondutors in all symmetry classes. Finally, in 
appendices we collect various exact Green's function identi- 
ties for superconductors, which are useful in the main text, 
and detailed derivations of several equations in the main text. 



II. CORRELATED CHIRAL TOPOLOGICAL 
SUPERCONDUCTORS IN TWO DIMENSIONS 

In this section, we will obtain topological order parameter 
for chiral superconductors. This topological order parame- 
ter is expressed in terms of zero frequency Green's function. 
This is analogous to the generalized first Chem number for 
quantum (anomalous) Hall insulators-'*, however, there are 
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differences between the two, among which the most promi- 
nent is the breaking of t/(l) gauge symmetry for supercon- 
ductors. More precisely, the Green's function matrix contains 
four submatrices, which are nonetheless not independent with 
each other The fermion operators are Qq., a - 1,2...N, 
where a refers to any degree of freedom other than mo- 
mentum. We can define the 2A^-component Nambu opera- 
tor as ^{k) - (ck,c[i^y. The Matsubara Green's function 

is defined as A;) = - e''^^Tr'i'a(T, /fc)*I'^(0, k) with 

a,/3 = 1,2... 2N. The Green's function are decomposed 
into four submatrices as 



g(iaj, k) 



GA(ioj,k) GsiicOyk) 
Gciioj, k) Go(ioJ, k) 



(1) 



where 



{GA)ap{iLO,k) = - J dTe"^\T,c,,{T,k)clp{Q,k)) 

{GB)al3{iC0,k) = - dTe""{TrCka{T,k)C-kl3{0,k)) 

Jo 

{Gc)ap{ico,k) = - dTe'"\T,c\jT,k)clp{Q,k)) 

(GDU(ico,k) = - f dTe'"'{Trc\^(T,k)c-kfi(Q,k)} (2) 
Jo 

Following Ref .-^^'^^ , we diagonalize the inverse Green's 
function as 



' (icij, k)\a{ia>, k)} — k)\a{iaj, k)) 



(3) 



We define those eigenvectors with fia{0,k) > as "R- 
zero", while those with fiaiO, k) < as "L-zero"^^-^^. From 
Eq. dAlOb . we know that the Green's function is Hermitian at 
zero frequency, more explicitly, we have 0~^(Q,k) - Q{Q,k), 
therefore, the eigenvalues of 0(0, k) is real. The spaces 
spanned by R-zeros and L-zeros are named as R-space and 
L-space. From the property of Hermitian matrix we can see 
that vectors within R-space are orthogonal to those within L- 
space, therefore, we can define Berry connections and cur- 
vature in the R-space. With these preparations, we can now 
propose one of the central results of this paper, namely a gen- 
eralized Chern number 



C, 



2nj 



d^kT, 



xy 



(4) 



where Tij - diJ[j - dj^,, and - -iYjR-space{k"\dk\ka), 
where \ka) are orthonormal vectors spanning the R-space. 
The simplest basis choice is \ka) - \a{i(jj - Q,k)). Eq.© 
strongly resembles the result for the quantum anomalous Hall 
insulators^, however, the R-zeros in superconductor cases 
are mixtures of particle/hole components. We also mention 
in passing that since only zero frequency Green's function is 
needed in Eq.®, we can also use real frequency Green's func- 
tion, which is the same as the Matsubara Green's function at 
zero frequency. Since we only need zero frequency Green's 
function, Eq.® can be evaluated efficiently by various ana- 
lytical and numerical methods. 



Next we will show the relation of Eq.(|4|i to another topolog- 
ical invariant defined for chiral superconductors, namely the 
winding number defined from interacting Green's function at 
all frequency, which is explicitly given as22id^ 

^'"i/ d^''dojTT[^''P(gd,g-')(gd,g-')(gdpg-')] (5) 

This topological quantum number measures the quantized 
thermal Hall effect^. We can prove the equivalence 



(6) 



by similar calculation to its analogue in the quantum Hall 
effect"*^. The key idea is to introduce a smooth deformation 
of giio), k) which does not change N2- Due to intrinsic prop- 
erties of Green's functions, a simple deformation connecting 
Q{ito, k) to an "effective noninteracting" Green's function can 
be found (see Appendix IB] for details) , it can thus be shown 
that N2 is fully determined by zero frequency Green's func- 
tion. This leads exactly to N2 - C\. The readers are referred 
to Appendix IB] for details of calculations, which is a straight- 
forward generalization of the derivations in Ref . . 

Similar to the case^'' of quantum Hall insulators, Eq.® 
shows that it is unnecessary to do frequency integral to obtain 
the topological invariant for superconductors, because zero 
frequency Green's function already contains sufficient infor- 
mation to determine the topological class. This greatly sim- 
plifies numerical and analytical calculations. 

The topological invariant given in Eq.lHli is equivalent to 
the Chem number calculated from an effective "noninteract- 
ing Bogoliubov-de Gennes Hamiltonian" 



h(k)^-g-\(jj^Q,k) 



(7) 



therefore, we may be tempted to think that Eq.® is a gen- 
eralization of the noninteracting topological invariants to the 
"renormalized energy bands", or "Bogoliubov-de Gennes 
quasiparticles". However, this picture is incorrect. Let us ex- 
plain more about this important point following Ref*^. In 
fact, in the "renormalized bands" or "quasiparticle" picture, 
we should use the self-consistent equation g^\a>a, k)\ua{k)} = 
0, which reveals the approximate quasiparticle poles. Taking 
advantage of the Dyson equation g^\cL>a,k) - u)a - hoik) - 
S(dLt„, k), we have 



[hoik) + Y(a)a,k)]\ujk)} = ojaMk)} 



(8) 



where ho(k) is the free part of the Hamiltonian and A;) 
is the self energy generated by electron-electron interaction. 
This "renormalized bands" or "quasiparticle" approach would 
suggest we define the topological invariant in terms of |m„(A:)), 
however, this approach is different from the zero frequency 
Green's function approach. Unlike this "renormalized band" 
picture, the zero frequency Green's function approach adopts 
a different equation'*'' 



[ho(k) + -L(0,k)]\ka} = ejka} 



(9) 



which follows from the Dyson equation '(0, k) - h()(k) + 
11(0, k). Eq.(|9]l is counterintuitive because the nonzero eigen- 
value e„ suggests we self-consistently use 1.{ea, k) instead of 
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2(0, k). As a comparison, we note that in Eq.®, \ua,(k)) has 
the clear physical meaning as the " quasiparticle"(or "would- 
be quasiparticle", because the concept of quasiparticle is most 
useful only when the lifetime is long), and the correspond- 
ing eigenvalues can be regarded as energy spectra. In con- 
trary to Eq.dS), in Eq.(|9) such clear physical meaning is ab- 
sent for \ka}, and the corresponding eigenvalues e„ cannot 
be interpreted as energy spectra, because of the lack of self- 
consistency in Eq.(|9l). 

To summarize the above discussions, the "renormalized 
bands" or the "(would-be) quasiparticle" approach prefers the 
self-consistent Eq.®, while the zero frequency Green's func- 
tion approach in this paper makes use of Eq.(|9). It has been 
emphasized'*^ that the self-consistent "renormalized bands" or 
"quasiparticle" approach following Eq.® is not suitable for 
calculating topological invariants, though it is an accurate tool 
for obtaining the energy spectra. On the other hand, the seem- 
ingly inconsistent zero frequency Green's function approach 
is an exact tool for topological invariants'*^, though it is a poor 
tool for the purpose of energy spectra. 



III. STRONGLY CORRELATED TOPOLOGICAL PHASES 
TRANSITIONS VIA GENERALIZED CHERN NUMBER 

Since our new formalism of topological order parameter 
can be applied to general interacting superconductors, we will 
use a simple ansatz of Green's function to illustrate its usage. 
Suppose that the Green's function of a 2d superconductor is 
given as 



Eq.® that 



g-\iu),k)^ 



lU) + n ■ T 

(w2 + n2)r 



(10) 



where n - (tlx, n,) - (sin k^, sin ky, m + 2- cos k^ - cos ky), 
and T - (Tv, Ty, Tj) are the Pauli matrices. In the — > limit, 
this ansatz becomes Q^^iio), k) x (iai + k ■ t + mT~)(aP' + k^ + 
rrp-y^ . If y = 0, this ansatz describes a mean field super- 
conductor with essentially free BdG quasi-particles. Gener- 
ally, this ansatz describes a strongly correlated superconduc- 
tor. When we replace mr. term by a k^j^ term, this Green's 
function becomes the "unparticle" propagator'*'*-'*'', with power 
law Green's function which is characteristic of conformal field 
theories. Unlike these gapless systems, our Green's function 
ansatz Eq.lfTOt describes a gapped system (for fermion), and 
the topological order parameter can be readily calculated us- 
ing Eq.®. 

Following the formalism in the previous section, we write 
down the zero frequency Green's function as 



(«2)r 



The Berry curvature T^- can be obtained as 



1 



(11) 



(12) 



where ff - n"/ + "2 "r follows from Eq.© and 



Ci(m) = 



0, m > 
-1, -2<m<0 

1, -4 < TO < -2 
0, TO < -4 



(13) 



It is very interesting to note that the highly nontrivial power 
y does not enter the topological number because only the 
Green's function eigenvectors are needed, which do not de- 
pend on y. The simplest way to obtain Eq.(fT3ll is outlined as 
follows. The generalized Chem number given by Eq.(|4) is un- 
changed as we tune to except at the three points to - 0, -2, -4, 
where the Green's function encounters singularities. Near 
these three singular points, we can linearly expand h", and ob- 
tain the changes of C]. For instance, ACi(m = 0) = Ci(to — > 
0^)-Ci(to — > ) = -1 can be obtained by linearly expansion 
n Si (kx, ky, to) near ^ = 0. 

The central point here is that the ansatz in Eq.lfTOli describes 
a strongly correlated superconductor, which can nevertheless 
be detected by the generalized Chern number for interacting 
superconductors in Eq.©. At the topological phase transition 
point, e.g. when to - 0, the situation becomes more interest- 
ing. Green's function becomes singular at this point, but it is 
does not fit into the non-interacting picture when y Q, since 
the small k behavior of Green's function is 



g(ia), k) 



{o? + k^y 

i<jj — k ■ T 



(14) 



which signifies a strongly correlated topological phase transi- 
tion. If 7 > 1 /2, Qiio), A:) — > when iiu), k) — > (0, 0), in sharp 
contrast to the non-interacting transitions with Qiiu, k) 00. 

It is worth mentioning that the same ansatz as Eq.lfTOll 
can be proposed for the quantum Hall effect, and an analo- 
gous strongly correlated topological phase transition can be 
obtained using the generalized TKNN invariant obtained in 
Ref,^ 



IV. CORRELATED TOPOLOGICAL 
SUPERCONDUCTORS IN 3D 

In this section we focus on time reversal invariant super- 
conductors in 3d, which in the non-interacting limit is classi- 
fied by an integer'"'**. In the Appendix lAl we show that the 
Green's function for superconductors at (/w, k) and (-/w, -k) 
is intrinsically correlated as illustrated by Eq. (|A7) . Time re- 
versal symmetry provides an additional symmetry of Green's 
function relating {iu, k) and (ioj, -k) as illustrated by Eq. (IC7l l. 
From these two equations we can obtain 



Y^{iu, kyL 



-Q{-ici), k) 



(15) 



where S = -iC^^T - Tv ® cr,,, in which t and <t are Pauli 
matrices in particle-hole and spin spaces [see the Appendix 
ICl . Since 2^ = 1, the eigenvalues of 2 is ±1, and we can 
choose a basis so that 2 is diagonal 



h 



-h 



(16) 
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Due to Eq.dTsTi, zero frequency Green's function in this basis 
is given as 



(17) 



The topological order order parameter of 3d time reversal in- 
variant topological superconductors is defined as the follow- 
ing winding number 



(even for noninteracting limit) for the angle from the pre- 
vious works. Therefore, the remaining problem is to evaluate 
in Eq.(fT9ll explicitly. A natural guess is that is proportional 
to W given in Eq.dTSll. namely = Wn. However, can only 
be defined mod In from the bulk information, because a In 
ambiguity can always be introduced depending on the details 
of interfaces between superconductors with different W. For 
instance, the simplest surface states of a topological supercon- 
ductors with W = 2 can be described by two surface Majorana 
cones, and the surface effective Hamiltonian reads ^^'^'-^' 



W = f d\€^"'''Ti:[(QdM-^)(QdyQ-^)(QdpQ-^)] (1^) „ V V r / , ,\ 



(20) 



= 1,2 k 



which has similar expression as its non-interacting analog 
however, Eq.JTSTl is well defined for superconductors in the 
presence of electron-electron interaction. As a byproduct, this 
shows that the integer classification of 3d time reversal in- 
variant topological superconductors is stable with respect to 
electron-electron interaction, although various possible frac- 
tional states are not included in this description. Therefore, 
free fermion superconductors in different integer class cannot 
be smoothly connected even if we add fermion-fermion inter- 
action effect. 

In the non-interacting limit, it is straightforward to check 
that Eq.JTSTl reduces to the winding number defined in terms 
of non-interacting Hamiltonianiii. 



V. GRAVITATIONAL TOPOLOGICAL FIELD THEORY 
AND GREEN'S FUNCTION AT ZERO FREQUENCY 

In this section, we focus on physical responses of 3d time 
reversal invariant topological superconductors and relate them 
to the Green's function. Since in superconductors the U(l) 
symmetry corresponding to fermion number conservation is 
spontaneously broken, we cannot get a topological field the- 
ory by coupling fermion to L'(l) gauge field, instead, we can 
get a gravitational topological field theory since the space- 
time metric couples to the energy-momentum tensor The 
gravitational topological term obtained in Ref.^^'^*' from chiral 
anomaly is given as 



rTrR A R 



lire 

224 j 8^ 



1536;r 



(19) 



where R"^^^ 



d^^vl^ + Ka^un - (/^ ^ and the diff^eren- 



tial form R = ^R^ydx'^ A dx^ . It is worth noting here that the 
pre-factor 1/2 in the first line of Eq.(fT9]l comes from the fact 
that the fermions are Majorana instead of complex. The form 
■^R A R is known as the Pontryagin class in mathematical 
literatures^*^. For general superconductors - 0{x, t) can be 
space-time dependent, but for time reversal invariant super- 
conductors, is quantized to be or ;r mod 2n. 

As we mentioned in SecU although the above gravitational 
topological field theory provides a description of interacting 
superconductors in this class, there is no explicit formula 



where rja^k - i^a,kT,T]a,ki)^ is a two component Majorana 
fermion operator satisfying 77^ — r]\j^. A mass term of the 
form Yia niarj^cTyria can induce a gap for the surface states. 
If mi,m2 > 0, the surface thermal Hall conductance^^ is 

K„ = 2 X "2^; on the other hand, if nii > 0, OT2 < 0, we 

jrk^T 

have /Cjy = X -5^. Therefore, the bulk topological invariant 
cannot fully determine the surface thermal Hall conductance. 
Following the relation^'''**^ between surface thermal Hall re- 
sponses and gravitational topological field theory, namely that 

6k~ T 

K„ - (assuming that the vacuum has - 0), we con- 
clude that - In for the case mi , m2 > 0, and = for the 
case nil > Q,ni2 < 0. To summarize the above calculations, 
if we know from the bulk that the topological superconductor 
has W = 2, we can at most conclude that = 2n mod 2jt, in 
other word, = is as equally possible as = 27r for this 
topological superconductor Therefore, we should modify our 
previous guess = Wn: to 



= Wk mod 27r 



(21) 



We will obtain this result by explicit calculation in Appendix 
|E] The Eq.(l2ni is among the central results of this paper 

Now let us find a different approach to calculate the gravita- 
tional angle. By analogy with the electromagnetic responses 
in topological insulators^, we propose the following expres- 
sion for the angle 



0^ dkodhj^^ dhTri^'p^'gd^.g-'gdyQ' 



xgdng-'gd^g-'gdrg-'] 



(22) 



where = ioj is the continuous Matsubara frequency( in the 
zero temperature limit), ki,k2, are spatial momenta, and 
is a dimensional extension parameter (Wess-Zumino-Witten 
parameter). The reference function g(ko,ki,k2,ky,k4 = n) is 
chosen as a trivial "flat-band" Green's function analogous to 
that in Ref.^'' [see Appendix |E] for explicit form of reference 
function for superconductors]. 

A direct Feynman diagram calculation of gravitational 
angle is much more involved than the electromagnetic an- 
gle, therefore we will justify Eq.(l22li as follows. In the non- 
interacting limit, Eq.(l22]l reduces to the Chem-Simons term 

,j{k)-^i[Ai(k),Aj(k)])-At(k)} 
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= ^ J d'k£'j''Tr[diAj+^iA!Aj]Ak (23) 

where A"^ = -i{if/'\dk,\ip^) are Berry connections defined 
in terms of Bloch states For a non-interacting Dirac 

fermion with Lagrangian 0'(/y''5^ + ot)^^'^^''*", from Eq.(|23]l it 
follows that 

e(m > 0) - 6(m < 0) ^ 1 ^^4) 
27r 2 

which is the same as the result obtained from chiral 
anomaly Since non-interacting Hamiltonian can al- 
ways be deformed to Dirac type at low energy, we conclude 
that Eq.(l22li is exactly the gravitational angle in the non- 
interacting limit. For general interacting superconductors, the 
quantization of 6 (due to time reversal symmetry) dictates that 
Eq.(l22b is the correct angle, which is analogous to the case 
of electromagnetic angle^'^^. The value in Eq.(l22li is of 
Z2 character because of the 2n ambiguity of the dimensional 
extension^^. This is consistent with the fact that the bulk of 
a superconductor can only determine the surface thermal Hall 
coefficient mod even integer^'', as we have discussed early in 
this section. 

Since we have justified Eq.(l22b as the correct gravitational 
angle, we can now derive Eq.dSTb from Eq.(l22li. The details 
of this calculation is left to AppendixlE] 

VI. GENERALIZATION TO GREEN'S FUNCTION 
SINGULARITIES IN GAPLESS SYSTEMS 

It is also straightforward to generalize the zero frequency 
Green's function topological invariants in this paper to gap- 
less systems, e.g. systems with Weyl points in momentum 
space. The zero frequency topological invariants are easier to 
use than the topological invariants with frequency integral^. 
For instance, a zero frequency Green's function topological 
invariant can be defined for a singular point (Fermi point, or 
Weyl point) in 3d gapless systems (insulator or superconduc- 
tor) as 




where we have chosen a small 2d sphere S around this sin- 
gular point, and S ' is its area element. The connection and 
curvature is defined as = -i TjR-spacAko!\dki\ka} and 'Fij = 
di^j - dj^i, where \ka} are orthonormal vectors spanning the 
R-space. Equation ( |25] ) is a direct generalization of Eq.(|4) to 
gapless systems. The construction of a 2d sphere around the 
singular point can be found in Ref."*^, in which topological 
invariants for Fermi points are defined in terms of frequency 
integral. Equation (|25] | has the advantage that no frequency 
integral is needed. For the definition of Eq.(l25l. it is crucial to 
realize that the zero frequency Green's function is a Hermitian 
matrix, so that a Chem number can be defined on the k space. 



VII. CONCLUSIONS AND DISCUSSIONS 

In this paper we have obtained several topological order pa- 
rameters for interacting superconductors. They are expressed 
in terms of zero frequency Green's function, and their equiva- 
lence to topological order parameters defined with frequency 
integral is explicitly shown. These new topological order pa- 
rameters are much easier to calculate in numerical and ana- 
lytical calculation. They produce topological invariants dif- 
ferent from the "quasiparticle" or "renormalized band" ap- 
proach, in which Green's function (or self energy ) at nonzero 
frequency should be self-consistently used. Taking advantage 
of the zero frequency Green's function approach, we analyze 
a nontrivial topological quantum phase transition beyond the 
non-interacting picture. The topological field theory coeffi- 
cient, namely the gravitational angle, has been expressed in 
terms of the zero frequency Green's function. 

For numerical calculation with finite size, at first sight 
the anomalous Green's function Gb,Gc would vanish be- 
cause of the fermion number mismatch. One viable ap- 
proach to circumvent this difficulty is to study the off- 
diagonal long-range order^^S. The idea can be outlined 
as follows. Since (0|cQ.(ri)c^(r2)|0) - (where |0) is 
the ground state) for any large but finite system, we can 
calculate (0|co,(ri)c8(r2)c|(r3)c],(r4)|0) instead. In the limit 
|ri - r2|, |r3 - r4| << |ri - rjl, this expectation value can 
be formally decomposed as (cQ.(ri)c/3(r2))(c|(r3)c^(r4)), from 
which (cQ,(ri)c8(r2)) can be extracted. In fact, this can be 
taken as the definition of {Ca(ri)c/j(v2)}- The only ambigu- 
ity is the global phase factor of Cooper pairing, which is not 
significant for our purpose. 

In the above simplified analysis we omit the t ( 
or w ) variable, and this variable can be straightfor- 
wardly added as (0|cQ.(ri,Ti)c/j(r2,T2)4(r3,T3)c|(r4,r4)|0) 
(or <0|cQ,(ri, icji)c0(r2, iW2)Cy(r3, !a)3)cj(r4, ia)4)\0) ), thus we 
can extract <Q,(ri, ri)c^(r2, T2)) (or <c„(ri, ;wi)c/3(r2, /W2)) ), 
whose Fourier transformation into the momentum space leads 
to the single particle Green's function Gb and Gc needed to 
define topological invariants. 

This is analogous to calculating magnetization in finite sys- 
tem. The naive expectation value of magnetization (in(r)) = 
0, but we can calculate {ma{ri)mt,(r2)} (a,b = x,y,z)- In the 
limit when |ri - r2| is large, we have the formal decompo- 
sition <OTn(ri)wi(r2)) ~ <mo(ri))<mfe(r2)), where {niair)) can 
be regarded as the definition of "magnetization" in large but 
finite systems. This definition is natural because {ma(r)) be- 
comes the true magnetization in the thermodynamical Umit. 
The nonzero magnetization in the thermodynamics remains 
visible as the long range correlations in large but finite sys- 
tems, though the latter do not have magnetization in a strict 
sense. 

In all the calculations in this paper, we have assumed that 
there is no ground state degeneracy (besides the phase of su- 
perconducting order parameter). For fractional phases which 
violate this condition, our formulas are not directly applica- 
ble. The extensions of the present approach to such exotic 
fractional phases will be left for future works. 
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We conclude with the remark that the zero frequency 
Green's function approach used in this paper is equally ap- 
plicable to topological insulators and superconductors in all 
symmetry classes in the "periodic table"" . The recipe is sum- 
marized as follows. We can obtain the zero frequency Green's 
function using analytical or numerical methods, then we can 
just calculate topological invariants using the "noninteract- 
ing Hamiltonian" h(k) = -G"'(w - 0,k) (for insulators) or 
h{k) - -0^^{ai - 0, k) (for superconductors). 

ZW would like to thank Xiao-Liang Qi and Xi Dai for dis- 
cussions. ZW is supported by Tsinghua University Initiative 
Scientific Research Program(No. 20121087986). SCZ is sup- 
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Appendix A: Exact Green's function identities for 
superconductors 

In this appendix, by straightforward calculation we will ob- 
tain several useful identities of Green's function for super- 
conductors without assuming any additional symmetry such 
as time reversal symmetry. Variant forms of some of these 
identities can be found in Ref.^^. The departure point of our 
calculation is the Lehmann spectral representation of the four 
submatrices in 0{ia), k) given as 



(GA) afi(icO, k) 

(GB) a(j(itO, k) 



{n\cka\m){m\cli^\ri} 



mn 
mn 



ioj - (E„ - E„) 
{n\cka\m){m\c-ki3\n) 
icD - (E,„ - E„) 

(n\c[i,Jm}{m\cl^\n) 



id) - (E,„ - E„) 

{n\c\)m){m\c-kii\n) 
ia> - (Em - E„) 



(Al) 



where \m) are exact eigenvectors of K - H - fiN with eigen- 
values Ei„ (H is the many-body Hamiltonian, fi is chemical 
potential, is particle number), D,„„ = e^^^ie^P^"- + e"^^"), 
and Q. is the thermodynamical potential defined by e"^" - 
jj-g-fiiH-fiN) _ -y^g ^jj calculations at finite temperature and 
then take the zero temperature limit, which is always implied 
in this paper. There is a important point hidden in the sim- 
plified notation in the above equations, namely that in the 
"anomalous propagators" Gb and Gc, the naive expressions 
such as {n\ciia\m}{m\c^iij3\n} vanish identically due to mismatch 
of fermion number, therefore, the precise meaning of the 
above formal Lehmann representations for Gb and Gc should 
be understood as follows. Let us focus on the zero temperature 
limit yS — > oo, then we have the naive expression 



{GB)aB(iOJ, k) 



z 

m 

z 



Q\cka\m){m\c-kp\'i) 
ioj - (E,„ - Eo) 

0\c-k/3\n}{n\cka\0} 
iu) + {En - Eo) 



(A2) 



where |0) is the ground state. This expression of Gb vanishes 
due to fermion number mismatch. The correct spectral repre- 
sentation should be modified as 



s /• ,x {0N-2\cka\m}{m\c-kfi\0N} 

(GB)a/i(lCi>,k) = 2j 

in 

-z 



i(o - {E,„ - Eq) 

{0N-2\C-k/3\n}{n\Cka\0N} 



iu) + (E„ - £0) 



(A3) 



where \Qf^) is the ground states with fermion number A^. Note 
that two ground states with fermions number and A^-2 enter 
this spectral representation. The ground state energy ( H-jiN) 
satisfies EoiN) ~ -£'o(A^ - 2) in the thermodynamical limit, 
which is simply denoted as Eq. Note that we have absorbed 
the chemical potential into the definition of Hamiltonian, in 
other word, E„ is the eigenvalue of H - fj.N. (If we define E„ 
as eigenvalues of H, we would have £'o(A^)-£'o(A^-2) ~ 2fi^'^). 
Eq.( IA3l ) should always be implied wherever the formal spec- 
tral representation given in Eq. dAH for Gb appears. Similarly 
modification of Gc is also understood in this paper 

With these preparations of Lehmann representation, we can 
obtain exact Green's function identities. The first exact iden- 
tity in this appendix is given as 



Go(/w, k) - -G^(-icj, -k) 
which can be obtained from the following calculation 



(A4) 



(GD)al3ii(^,k) = VA 



mn 



mn 



{n\c_i^Jm}{m\c-ki3\n} 
{m\c-k/}\n}{n\c[jm) 

ICO t^mn 

{m\c-kp\n){n\c[i^J,m) 



= -(GA)fia(-i<^,-k) 

where we have defined E„„ = E„, - E„. We mention that an 
variant form of Eq. (IA4) has been given in Ref.^^. 
The second identity is 



GB(i(^,k) - -G]j(-ia>, -k) 
which can be calculated as follows 

^ rx {n\cka\m){m\c-ki3\n) 

(GB)ap{l<^,k) = > D,nn- 



(A5) 



mn 



mn 



ico E/jjij 
{m\c-ki3\n}{n\cka\m} 



-ioj - E„, 
(GBha(-i(^, -k) 



We can check that Eq. dASI l can also be obtained directly from 
the more precise spectral representation given in Eq. (IA3) . 
By analogous calculation we can also obtain 



Gciito, k) — —G^(—ia>, —k) 



(A6) 
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Now we define a charge conjugation matrix 

1 



1 



then we can readily obtain that 

Cg(ito, k)C-^ = -g^i-ico, -k) (A7) 

which can be calculated as 

nnr i\n-i ( GD{i(^,k) Gc{ioJ,k) 
CG{ioj,k)C = [G^(ioj,k) G^ick) 

_ I -G^{-ia>, -k) -G^{-ia>, —k) 
\ -G^(-ia>, -k) -G^(-ia>, -k) 

= -g^(-iaj,-k) 
By similar calculations, we can also obtain that 

GB(icj,k) ^Gl(-iLo,k) (A8) 

and 

GA{iciJ,k) =Gl{-iuj,k) (A9) 
from which it follows that 

g{ia),k) ^ g\-ioj,k) (MO) 

Appendix B: Derivation of A'2 = Ci for superconductors 

The key idea of this derivation is analogous to the one given 
in Ref.^^, though the calculations for superconductors is a lit- 
tle more involved because of appearance of four submatrices 
Ga,Gb,Gc,Gd in the Green's function Q{i(jj,k), as we will 
show below. 

The key idea is to introduce a deformation of Q which 
smoothly connects Qiiu, k) and an "effective non-interacting" 
Green's function (/w + ^(0, fc)"')"' as 

g(iLo, k,A) = (l- A)g(ico, k) + Aiioj + g-^O, k))-^ (Bl) 

where A e [0, 1]. It is not evident that this deformation is 
smooth, since in principle ^(/w, k. A) can have zero eigenval- 
ues. It turns out that zero eigenvalue does not appear because 
of intrinsic properties of Green's function as we show below. 
It is convenient for our purpose to decompose = 01+ 102, 
where 0u02 are both Hermitian matrices. The explicit form 
of 02 is given by its four submatrices 

(G2A)a/i(ico, ^) = - ^ d„,„{n\cka\m){m\cl^\n) 

inn 

(G2b)q-/j('w, ^) = - ^ d„,n{n\cka\m){m\c-k/i\n} 

mn 

(G2c)a/}(i<^, ^ ~ X '^'""^"I'^-iQ^I'")^'"!'^*/?!") 

mn 

{G2D)aii{i(^,k) = -^d,n„{n\c\j^m){m\c-kp\n) (B2) 



where c/,„„ = oj'he' -e y- ' defined in the Ap- 

pendixlA] The crucial property of 02 is the following. Taking 
an arbitrary 2N component vector = 1^'))"'^, in which 
\b) are A^-component vectors, we have 

<(^l^2l'/'> = - ^ d,„„[\Am„\^ + A*„„B,„„ + A,nnBl,„ + |B„„/] 

nm 

— ^ ^ ^ dmn\Amn + Bmnl (1^3) 
mn 

Here we have defined Am„ = 2q. fli^(n|c<:Q.|m) and B,„„ = 
2a bain\c^^i^Jm} (here Oq,, ba are the components of \a), \b}). 
The following calculation is almost identical to the case of 
insulators^*, which we will briefly outline. From Eq. (IB3l l it 
follows that 



sigR{{t//\02{ico, k)\iff}) = -signw 



(B4) 



for an arbitrary \ifr}. Suppose that is an eigenvector 
of 0(iLj,k,A) with eigenvalue ju^\ then it can be obtained 
that \m{n-^\ioj,k,m = <tA|tA)-'[(l - A){4,\02{iojM<lj) - 
AcijYjsVI^ s?'(Mi^ + f?) where we have expanded = 
Yjs^s{i<J^,k,A)\s{k)), with \s{k)) being the eigenvector of 
' (0, k) with eigenvalue e,. From this equation and Eq.( IB4l i 
we can see that when icj + 0, Im[/i^'(/w, fc, /I)] + 0. On the 
other hand, when iuj = 0, ^(0, k, /I) = ^(0, k) is independent 
of A, and Re[yU^'(/w, k, X)\ + 0. Summarizing the above calcu- 
lations, we have pT^iiu, k. A) + 0. Therefore, the deformation 
given in Eq. dBll l is smooth, and we have N2{A = 0) = N2{A = 
1). Now we can calculated N2 = N2{A = 0) by calculating 
N2{A =1) with "effective non-interacting" Green's function 
{id) -H ^(0, fc)"')"'. By a straightforward calculation, we are 
led to Eq.©. 

The interested readers are referred to Ref.— for a similar 
derivation for the case of insulators. 



Appendix C: Exact Green's function identities in the presence of 
time reversal symmetry 

The basic calculational tool in this section is the Lehmann 
representation and the following identity 



{a\b) = {-\f'{b\a), {a\b) = (-lf"{b\a) 



(CI) 



which follows from t^\a) = (-l)'*'"|fl), where is the 
fermion number of \a), and |a) = T\a). It is worth noting that 
Eq. dCll l are identities for exact many-body eigenstates \a} and 
1^), which should not be confused with the analogous identit- 
ties for non-interacting Bloch states, in which (-1)'*' factor is 
replaced by - 1 . Taking advantage of time reversal symmetry, 
we will obtain several exact identities. The first one is 



Gliicj,k) = TGA(ico,-k)T'^' 



(C2) 



where T = /cr,, is the time reversal symmetry matrix defined 
by r = icTyK {K is the complex conjugation operator). The 
transformation of fermion operator is TckaT^^ = Tj/iTapc^kp- 
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Eq.( IC2l ) is analogous to the time reversal identity in time re- 
versal invariant insulators^^'^^. Eq.( |C2| ) can be derived explic- 
itly as 



{n\ci,a\m){m\cl^\n) 

mn ""' 



mn 



{n\cka\m){m\c[ij\n) 

ICO I^iiw 



{n\c_i,s\m}{m\c''_i^^\n) 

= }_^t:^tp,d„,„ -— ^- 

mn ""' 

= (TGA{iaj,-k)T'''')pa 

In the third line, we have used the first one in Eq. dCH 
by taking \a} = |n) and \b) - Qq.|ot), which leads to 
\b) = TciiaT^^T\m) - TayC_ky\m), and finally <n|cfa,|m) - 
{-\)^'"^^Tly{m\c'^ ^,^\n). By a similar calculation, we can also 
obtain 

Go(ia), k) = T'GDiioJ, -k)T (C3) 
The third identity is 

Gliioj, k) = -T^'GcdtJ, -k)T'^ (C4) 
which can be explicitly calculated from 



k) = y,A 



mn 



mn 



mn 



{n\cka\m){m\c-ki3\n) 
lU) - h,„„ 

{h\cka\m}{m\c-ki3\n} 

10) t^mn 

KrTps{m\c[jn}{n\cl^\m} 



lOJ — i^inn 

Ky^6'^n\c^J^){m\c\^\n) 



10) E„„ 
^ -Tls^Gc)6yii(^,-k)T'l^ 
= -{T^Gc{io),-k)T^)pa 

in which we have used Eq. dCll ) several times and the identity 
ji _ which follows from r'T = 1 and T*T - -1. 

Similarly we can obtain 



G^(io),k) — -TGB(io),-k)T 



(C5) 



Now we can define a time reversal matrix in the Nambu space 

as 



T = 



(C6) 



Tg(io),k)T- 



With this definition, we can calculate 

TGA{io),k)T^ -TGB{io),k)T 
-r"'"Gc(/w,fc)r"'" T'GD(ioJ,k)T 

G^(io), -k) G^(io), —k) 
G^(io), -k) G^(io), -k) 

= g''(io),-k) 



which is written compactly as 

rg(io),k)T-^ = g^{io),-k) 



(C7) 



Appendix D: Green's function in the Majorana operator basis 

Since the Majorana basis is convenient for some problems, 
we will present several identities in this basis for references. 
Let us start with notations. The fermion operators are referred 
to as Cka, where k is spatial momentum, a - 1,2 .. .N refers 
to any other degrees of freedom including spin, orbital, etc. 
We also define 2N Majorana operators j^^, /3 - 1,2... 2N by 

Jka = Qff + ci^^, 7k,N+a = KCka ' c'_^„)\ (ff = 1, 2 . . . A^XDl) 

Since our results do not depend on the basis we use, we will 
formulate them in the most convenient basis, namely the Ma- 
jorana basis defined above. The Matsubara Green's function 
in this basis is defined as 



Qapij, k) = -{Tryka{T)ykl}{Q)) 



(D2) 



where is the imaginary time ordering. The fre- 
quency domain Green's function is defined as G(io)„,k) = 
Jq dTexp(io)„T)G(T,k). The Lehmann representation of 
Green's function at finite temperature is given by 



, {n\yka\m){m\y-ki3\n) 
io) - (E,n - E„) 



-HE,, 



)(D3) 



where p is the inverse temperature, \m) are exact eigenvectors 
of K - H — fxN (fi is the chemical potential, = Yjka '^la'-ka 
is the fermion number operator), and Q is the thermodynamic 
potential defined by e"^" = Xaking advantage of 

this spectral representation, we can obtain several exact iden- 
tities. The first one is 



0'{io),k) = 0{-io),k) 



(D4) 



which is the same as the identity in insulators. The second one 
is 



0{io), k) - -Q^{-io), -k) 



(D5) 



which has no counterpart in insulators. This identity signifies 
major difference between the superconductors and insulators, 
and will be crucial for construction of interacting topological 
order parameters for superconductors in 3d. It can be appreci- 
ated that in contrast to the case in insulator, superconductors's 
Green's function at k and -k are intrinsically related because 
of the cooper pairing between opposite momenta. 

In the zero temperature limit, the discrete Matsubara fre- 
quency variable o) becomes continuous, and the Lehmann rep- 
resentation of the Matsubara Green's function reads 

s Y^r^O|yio-l»j)<»j|r-t/s|0) <0|r-i/j|m)<yt„|0> 
^ 10)- (E„, - Eo) 10) + (E,„ - Eo) 

where |0) is the ground state. 

We mention in passing that the derivations in Sec|B]can be 
simplified if we use the Majorana basis, in which the compli- 
cation of four submatrices in Eq. (IB2t can be avoided. 



9 



Appendix E: Gravitational 9 coefficient and zero frequency 
Green's function 

In this section we derive Eq. (l2Tl i. which relates the grav- 
itational 6 coefficient to Eq.lfTSll. First, we extrapolate the 
Green's function ^(A:o = icL>,ki,k2,ki) to 0(ko,ki,k2,ki,k4), 
with the reference Green's function at A;4 - n chosen as 
GA{ia),ki,k2,ki,n) = -j^lAfxAf, GD(ioJ,ki,k2,k2,7T) = ^Inxn, 
GB(i(ii,ki,k2,kT,,7T) = jT'^', GciicL>,ki,k2,k3,n) = jT, where 
A is a positive number with dimension of energy, and T = /cr,. 
is the time reversal operator in the spin space. Physically this 
Green's function describes a mean field trivial s-wave super- 
conductor with momentum independent pairing. 

By an calculation analogous to that given in Ref.^^ [see also 
Appendix iBll. we can deform 0(i(i),k) to (io) + 0^^{O,k)Y^ 
without encountering singularities. This fact will be crucial 
for the entire calculation, because it implies that zero fre- 
quency Green's function can fully determine the topological 
invariants. 

Taking advantage of the time reversal symmetry Eq. (IC7) . 
we can also write 6 in Eq.(l22]i as 

xgdpg-'gd^g-'gdrg-'] (ei) 

in which g in the k^ e [-tt, 0] region is determined 
by Tg(iaj,kuk2,k3,k4)T-^ = ^^(/w, -^i, -^2, -^3, -^4) 
and Cg(ia>,ki,k2,k3,k4)C^^ = -0^ {-ito,-k\,-k2,-k-i,-k4). 
Eq. dEU is evidently unchanged in smooth deformations. Af- 
ter deforming the Green's function to (/w + g^^{Q,k)Y^ in a 
similar fashion to Appendix iBl by a direct calculation we can 
express 6 in Eq.(l22b in terms of zero frequency Green's func- 
tion explicitly as 

OI2n = -L^J^^dk4j cPke'^^tATiiTki] 



= CS{k4 = 0) - CS{k4 = n) 
= CS(fe4 = 0) 



(E2) 



in which the Chern-Simons term 



CS = J d\e'''Tx{{Tij{k)-h[mi{k),mj{k)])-mk{k)] 

" i/ f''''Tr[5/^; + ^'^/^,]^* (E3) 



where 



mfik)^-i{ka\ — \k/3) 
oki 



in which \ka) is an orthonormal basis of the R-space spanned 
by R-zeros [see the main text, below Eq.lO, or see Ref.i^] 
. In the k integral in CS term, k4 is fixed to be constant, and 
d^k - dkidk2dk-i. We have used the fact that CS(A;4 = ;r) = 
because g{k4 = ;7r) is a trivial reference function. 



Now we can proceed to calculate the Chem-Simons term, 
which is at this stage already expressed in terms of zero fre- 
quency Green's function. We will show that we are led exactly 
to Eq.dTTTi by direct calculation. The physical Green's func- 
tion at zero frequency has the following form as obtained in 
the main text 



Q{kuk2,kj,) 



(E4) 



or 



.-(0,.,*.,. = )(E5, 

in which ^4 = is implied. From the chiral symmetry 
l.g{Q,k)Yr^ - -g(0,k) it follows that the eigenvectors \a} - 
(u, v)^ and \/3} = (u, -v)^ of g~^(Q, k) form chiral paks as 

g-\o,k)\/5) = 



with 



All R-zeros (namely eigenvectors \a(0, k)} with fia(0, k) > 0) 
span the "R-space"^'' at each k, on which the Chem-Simons 
term is defined. Due to the specific form of we can find 
all the R-zeros by finding the eigenvector of Q^Q. Suppose 
that iJ^Q'''Q\v} = |y) with <v|v> = 1, then la) = (±fiQ\v), \v)/ is 
an eigenvector of g^^ with eigenvalues iju. Therefore, a basis 
of R-zeros has the following form 



\a(0,k)}. 



^|2 \ l^«> 



(E6) 



with yUfj, > 0, where \va) is an component column vector 
satisfying yu^<3^(3|v„) = \va) and normalized by (vq.|v„) - 1 . 
The magnitude of eigenvalues //q, is irrelevant for calculations 
of topological numbers, so we deform all //q, = m > 0. As a 
result, the eigenvectors \a{Q, k)) in Eq. (IE6b are deformed to 



(E7) 



where is a dimensionless matrix satisfying H'li - 1 because 
/i^,<3'(3|v„) = |vv) is preserved during the deformation. To cal- 
culate the Chem-Simons term over the R-space in Eq. (IE3) , we 
can choose a trivial basis in which |v„) is A:-independent (this 
is possible because the fiber bundle spanned by \va) is trivial) 
, and the R-space Berry connection in this basis is given as 



(E8) 



where the crucial 1/2 factor comes from the 1 / V2 in Eq. dETI ). 
otherwise y[, would become a pure gauge. Eq. dESb is not a 
pure gauge just because of the crucial 1 /2 factor. This is the 
key of calculations in this section. 
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It is worth noting that although Chem-Simons term is basis 
dependent, a basis change (a "gauge transformation") can at 
most change the Chern-Simons term by an integer Our central 
result Eq.dSTb is a mod In equation, which is thus not affected 
by basis choices. The Chern-Simons term in Eq. (IE3b can be 
calculated directly from Eq. dESt as 

CS(fe4 = 0) = — ^ r d\ff''PTr[{'R-^dJ<){'Rr^dy'R){'R-^dp'R)] 

4871^ J 

= —Win) (E9) 



which is equal to -jWiQ) because smooth deformation can- 
not change the winding number. The minus sign comes from 
the exchange of H and compared to the definition in 
Eq.lfTSll. Summarizing the above calculations, we can see that 



— = CS(;fc4 = 0) = - WiQ) (mod integer) (ElO) 
2n 2 

which is exactly the Eq.(|2TJ. 
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